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Abstract 

Neutron electric dipole moment (EDM) due to single quark EDM and 
to the transition EDM is calculated in the minimal supersymmetric stan- 
dard model. Assuming that the Cabibbo-Kobayashi-Maskawa matrix at 
the grand unification scale is the only source of CP violation, complex 
phases are induced in parameters of soft supersymmetry breaking at low 
energies. Chargino one-loop diagram is found to give the dominant con- 
tribution of the order of 10 -27 ~ 10~ 29 e-cm for quark EDM, assuming 
the light chargino mass and the universal scalar mass to be 50 GeV and 
100 GeV, respectively. Therefore the neutron EDM in this class of model 
is difficult to measure experimentally. Gluino one-loop diagram also con- 
tributes due to the flavor changing gluino coupling. The transition EDM 
is found to give dominant contributions for certain parameter regions. 
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1. Introduction 



The supersymmetric theories now stand as the most promising candidate for the 
unified theory beyond the standard model []]]. The supersymmetry helps to resolve 
the gauge hierarchy problem [0. Moreover, the accurate data favor remarkably the 
supersymmetric grand unified theory (GUT) over the nonsupersymmetric theory 
||]. This fact has fulfilled the promise that accurate measurements of coupling 
constant strengths at low energies can distinguish various alternative candidates for 
the grand unified theories by extrapolating the renormalization group trajectories 
to higher energies. 

Among many problems in particle physics, the violation of CP invariance is one 
of the phenomena that are least understood. The primary reason for this unsatis- 
factory situation is that the experimental verification of the CP violation is so far 
limited to neutral Kaon dacays into two pions. We expect to obtain more experi- 
mental informations on the CP violation from the B-factory soon. The CP violation 
is not only important as a fundamental symmetry property, but also needed to ex- 
plain the cosmological baryon asymmetry of our universe ||. Therefore it is most 
desirable to have additional experimental informations on the CP violation. Apart 
from the forthcoming experiment with the B-factory, we have one more promising 
observable for the CP violation: the electric dipole moment (EDM), in particu- 
lar those of neutron and electron || . Since we can hope for further improvements 
of experimental precision, especially for that of neutron, we expect that the EDM 
will provide a precious clue of the CP violation. The CP violation in the minimal 
standard model arises solely from the Cabibbo-Kobayashi-Maskawa matrix of the 
Yukawa coupling constants of the Higgs field. Therefore progress in the study of 
the CP violation provides important informations on the Higgs field which is most 
elusive in the standard model. 

In supersymmetric models, we have more possibilities for complex parameters 
beside the Cabibbo-Kobayashi-Maskawa matrix, even if we have only the minimal 
particle content of the supersymmetric standard model. These complex parameters 
become additional sources of the CP violation. Among the parameters of the su- 
persymmetric models, those associated with the soft breaking of supersymmetry are 
least understood. The early studies of the EDM in supersymmetric models have 
revealed that generic complex parameters for the soft breaking give too large EDM 
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unless superpartners in the loop are heavier than at least a few TeV ||. Although 
this type of models are not excluded, it is perhaps more attractive and natural if we 
can control the phases of the soft breaking parameters so that superpartner masses 
of the order of the electroweak scale are naturally allowed. There have been many 
studies on this issue |7|, ||. 

In the most popular model, the supergravity with the hidden sector provides a 
definite pattern of the soft breaking of supersymmetry [jl], [|. If we assume a simple 
model for the hidden sector and the supergravity couplings, we obtain that these 
parameters are real at the grand unification scale or the Planck scale. However, 
the soft breaking parameters that are really manifest at low energies will become 
complex since nonvanishing phases will be induced by the renormalization group 
flow involving the Cabibbo-Kobayashi-Maskawa matrix |L(], |!|. Moreover, flavor- 
changing couplings of gluino are also induced [IT|. This radiative effect becomes 
important when the Yukawa couplings are large. It is worth examining the CP 
violation due to this radiatively induced phases of the soft breaking parameters, 
since it is now certain that the top quark is quite heavy |[2| and requires a large 
Yukawa coupling. 

More recently, there have been a number of studies on the neutron EDM in the 



supersymmetric models El or in two Higgs doublet models JT4 



In the nonsupersymmetric minimal standard model, it has been proposed that 
the transition quark EDM can be more important than the single quark EDM to 
explain the neutron EDM ]TJ, [R| [T7|]. 

The purpose of this paper is to examine the neutron electric dipole moment in 
the minimal supersymmetric standard model. We assume that there are complex 
parameters only in the Cabibbo-Kobayashi-Maskawa matrix at the grand unification 
scale. We examine the effect of the radiatively induced phases of the soft supersym- 
metry breaking parameters on the neutron EDM. In performing the renormalization 
group analysis, we have taken account of the effect of gaugino masses together with 
the universal scalar masses. We shall also consider the transition quark EDM in 
the supersymmetric models. We find that the single quark EDM is of the order of 
10 -27 ~ 10~ 29 e-cm for the light chargino mass to be 50 GeV. We also find that 
the transition EDM is of the order of 10 -25 ~ 10~ 27 e-cm and hence contributes the 
neutron EDM of the order of 10~ 29 ~ 10~ 31 e-cm, if we take account of the large 
QCD enhancement due to the penguin diagrams. In both cases, the neutron EDM 
in this class of models is too small to be detected in forthcoming experiments. 
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In sect. 2, we introduce soft breaking parameters of supersymmetry and scalar 
particle mass matrices. In sect. 3, we analyze the single quark EDM. There are 
two classes of contributions: the chargino loop and the gluino loop. In sect. 4, we 
examine the transition EDM. Appendix is devoted to describe our results of the 
renormalization group equations and our inputs. 

2. Soft breaking parameters of supersymmetry 



We consider the minimal supersymmetric standard model (MSSM) which con- 
tains left-chiral supermultiplets for three generations of quarks (U c , D c , Q) and 
leptons (E c , L), gauge bosons of the SU(3) c xSU(2)LxU(l)y and two Higgs dou- 
blets (H u , H d ). Boldface letters such as Q denote vectors in generation indices, 
and the suffix c denotes the antiparticle. The supergravity with the hidden sector 
provides the pattern of the soft breaking of supersymmetry [p], ||. In simple models 
of this type, the soft breaking parameters are real and universal at the grand uni- 
fication (GUT) scale or the Planck scale, and a complex phase appears only in the 
Cabibbo-Kobayashi-Maskawa matrix at the scale. Renormalization group running 
induces complex phases in the soft breaking parameters at the lower scale. 

One can write superpotential of MSSM as follows, 

-W = U cT Y u H u ■ Q - D cT Y d H d ■ Q - E cT Y e H d ■ L + fiH u ■ H d , (2.1) 

where boldface letters Y are Yukawa couplings as a matrix in generation indices. 
The inner product of SU(2) indices is abbreviated by the • as H u -H d = {ia2H u ) T H d . 

Soft supersymmetry breaking is given by the following terms in the Lagrangian: 
1. scalar mass terms, 

-EmJ^ (2.2) 

i,3 

The 0j's are all the scalar particles and m|„- is a hermitian matrix. In the 
supergravity-induced models, this will be universal at the GUT scale, i.e., 
m?- = m 2 5ij at the GUT scale. It runs by renormalization group flow. 
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2. A terms (trilinear scalar couplings), 



u R A u h u ■ q L + d R A d h d ■ q L + e R A e h d ■ £ L + H.c. (2.3) 



In the supergravity-induced models, A's are proportional to Yukawa couplings 
at the GUT scale, A = AY. The A terms run by renormalization group flow. 

3. B term, 



4. gaugino mass terms, 



Bfjihu ■ h d + H.c. (2.4) 



-^MA^ + H.c. (2.5) 

i=i 



If GUT is embedded in the supergravity-induced models, Mj (i = 1,2,3) 
are universal, i.e., Mj = M at the GUT scale. Gaugino masses Mj run by 
renormalization group flow. 

We can write scalar-quark mass terms in the following form: 

- (*l *« M l (Z) - (al al) «S (t) • (2 - 6) 



where wl,/? (d>L,R) is a u- type-scalar-quark (d-type-scalar-quark) field column vector 
in generation indices. 

Denoting the scalar-quark mass matrices at the GUT scale with the suffix 0, one 
can find 

2 (ml L l + M{M u (A + pLC0tp)M[\ 

"° " ^M u (A + //cot/3) m| R l + M U M.\ J ' l ' J 

2 (ml L l + M\M d (A + / ,Un(3)M d \ 

d0 " ^M d (A + //tan^) m\ R l + M d M\)> { ' } 

where M u and M d are matrices in generation indices for u-type and d-type quark 
masses respectively and are given by Yukawa couplings Y and vacuum expectation 
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values v of Higgs fields, 

M u = Y u ^=, M d = Y d ^=. (2.9) 

The tan/5 is defined as v u /v d . The generation independent part of the scalar-quark 

masses are given by the universal mass and the contribution from the D-term [|T^] : 

m 2 L = m 2 -M/(-cos2/3)Q-jjsin 2 ^), (2.10) 
2 

m \ R = m 2 - -M z 2 (-cos2j3)sm 2 6 w , (2.11) 
3 

m\ L = m 2 + M z 2 (-cos2/5)Q-isin 2 ^, (2.12) 

1 | B = m 2 + ^M z 2 {-cos2(3)sm 2 6 w . (2.13) 



m 

c 

In writing the above formulas, we assumed for simplicity a universal form of 
the supersymmetry breaking, namely the universal scalar mass m and the trilinear 
scalar coupling with the universal parameter A. Moreover, all the parameters are 
real except the Yukawa coupling constants which appear in the quark mass matrices 
at the GUT scale. Therefore the Yukawa couplings or the Cabibbo-Kobayashi- 
Maskawa matrix is the only source of CP violation in this model at the GUT scale. 

The renormalization group flow changes these scalar-quark mass matrices at the 
lower scale [jn|. At the lower scale, these mass matrices can be written as follows: 

2 _ (ml L l + MiM u + Sm 2 L (A u + /icot/31)Mj \ ( . 
• ^ ! M u (At + Aicot/91) m\ R \ + M U M\ + Sml R } '■ ~ 1 



d \ M d (A d + fi tan (31) m 2 - fi l + M d M\ + <5m 2 - R y ' l ' ; 

where <5m 2 's are hermitian matrices defined from the soft SUSY breaking scalar 
mass parameters which are determined from the renormalization group flow and 
have off-diagonal elements at the electroweak scale. After the renormalization group 
running, the parameter A's become matrices in generation indices and the left- 
left (LL) and right-right (RR) blocks have off-diagonal elements. The off-diagonal 
terms of the matrices A's turn out to be important for the neutron EDM. We show 
the renormalization group equations and their solutions for the matrices A's in 
Appendix [A]. 

It is convenient to rotate the scalar-quark wave functions by the same amount as 
to diagonalize the mass matrices of the quarks themselves, although scalar-quarks 



2 l + M ] d M d + 8m\ L (A d + /itan/?l)M t 
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are not in mass eigenstates in this basis. This basis of the scalar-quark wave function 
is usually called super KM basis. Denoting the wave functions in the super KM basis 
with a prime, we obtain quarks and scalar-quarks as 

(2.16) 
(2.17) 
(2.18) 
(2.19) 

where the quark mass matrices are diagonalized in generation indices with the ma- 
trices V u l and so on: 

V uR M u V{ L = m u , (2.20) 
V dR M d V\ L = m d . (2.21) 

The Cabibbo-Kobayashi-Maskawa Matrix is defined as K = V u l V\ l . We use 
a parametrization and explicit values of K as given in Eqs. ( |A.35| ) - ( |A.37 ) of 
Appendix [Sj. We also rotate 6m 2, s and A's as follows: 

.2/ 




5mf H = V fH 6mj H V} H , (2.22) 
i' f = V fL A f V\ L , 



A> = V fL A f V\ L , (2.23) 



where f — u,d and H = L,R. 

Scalar-quark mass matrices in this basis are given by 



M' 2 




ml T 1 + ml + 8ml 



'"uL* ~ "*u ~ " " L uL 1 r>~y» r*-J ,, "u j (2 24) 

m u (Al + ficotpi) ml R l + ml + 8ml R r 



M? - ' VdL ° 



d ~ \ V dR d \ 




'm| L l + m\ + Sm^ L (A d + //tan/31)ra d ' 
n m d {A d + titan (31) m\ R l + m d + Sm~ R/ 



(2.25) 
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Here we neglected the primes in the right hand sides. Hereafter we consider in this 
basis. We show A's at the electroweak scale in this basis which are solutions of 
renormalization group equations in Appendix [A|. We give initial conditions for A's 
at the GUT scale Mqut as in Eq. ( |A.28|) . The A's are diagonal at the GUT scale. 
Then A e is diagonal at all the scale, but A u and Ad are not diagonal at the lower 
scale. It is convenient to separate A's into two parts, 



A L/ + A 



Mf, 



(2.26) 



where / = u,d,e. The second term Am/ is proportional to the universal gaugino 



mass M and satisfies the same renormalization group equations ( |A.22 ) - ( A. 24 ) as 
Af. Their initial conditions at the GUT scale are 



A Mm (0) = A Md {0) = A Me (0) = 0. 



(2.27) 



The first term A^f satisfies linear equations (|A.31| ) - ( |A.33|) which are obtained by 
deleting gaugino masses in Eqs. ( |A.22| ) - ( |A.24| ). Our results of the renormalization 
group analysis are summarized in the Appendix by giving the matrix Af at low 
energies in Eqs. flQgp - ([031). 



3. Single quark electric dipole moment 



The matrix element of the electromagnetic current j' M (0) between a Dirac fermion 
with momentum p and spin component s can be written in terms of four independent 
form factors F; as follows: 



(Pf,Sf\j^(0)\p t ,Si} = u(p f ,s 



f) 



^F 1 (q 2 )+ia fiU q l 



2m 



+ fe 7 



2m 

u(Pi,Si), (3.1) 



where q = pj — Pi- The electric dipole moment (EDM) d of the spin-| particle is 
given in terms of the form factor F^{q 2 ) as 



(3.2) 



since a Dirac particle with spin vector s and with the EDM d interacts with the 
weak external electric field E as 



J 'mt 



2ds ■ E. 



(3.3) 



This interaction violates the CP invariance. 

We calculate the neutron EDM in the MSSM. Contrary to the minimal standard 
model without supersymmetry, we have contributions already at one-loop. There 
are two classes of contributions: 

(a) chargino loop, as shown in Fig. |l] 

(b) gluino loop, as shown in Fig. |2| 

Although the complex phase is assumed to be only in the Cabibbo-Kobayashi- 
Maskawa matrix at the GUT scale, the matrices Sm 2 L and A in the scalar-quark 
mass matrices also have complex phases at low energies as discussed in the previous 
section. We evaluate the contribution of these radiatively induced phases to the 
EDM. 

There are a few diagrams which give significant contributions to EDM. Especially, 
the diagram of Fig. [I] gives a dominant contribution for EDM since it involves the 
Yukawa coupling constant of the top quark most directly. The EDM from the 
diagram in Fig. |T| can be written as 



cL = 



e 2 

— — : -(tan/3 + cot/3)— m d Im 
v z 



KUF 



M„ 



LR m * ) hh 



(3.4) 
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where the subscript LR means the left-right block of the scalar-quark mass matrix 
(Eq. (|2.24j )) and the subscript 11 denotes the (1,1) component in the generation 
indices. The subscript hh denotes the (2,2) component in the wino-Higgsino mass 
matrix M x which is given by 



-(w, 



f! 



K R )M X 



w L 
h 



dL 



H.c. 



(3.5) 



Mo 



V2M W sin 



\/2M w cos f3" 



(3.6) 
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The function F(x) is given by 



F(x) = 



6(1 -x) 



-(5 - 12s + 7x 2 + 2x{2 - 3x) logx). 



(3.7) 



To evaluate the expression Eq. (|3.4j ), we need to diagonalize the scalar-quark mass 
matrix M~ 2 and the wino-Higgsino one M x . The former is explicitly diagonalized by 
a 6 x 6 unitary matrix Ua, 

UuM'M = Ml (3.8) 

The Mj is a diagonal matrix whose elements are all positive. On the other hand, 
the latter can be analytically diagonalized by two orthogonal matrices Ur and Ui 
as 



U R M X U T L 



m x i 
m x 2 



where the mass of light (heavy) chargino is denoted as m x i (m x 2 ) 
m x i, m x 2 = - J{M 2 - /i) 2 + 2M$y(l - sin 2(3) 



T V /(M 2 + /i) 2 + 2M 2 ,(l + sin2/?) 
Using those matrices, Eq. (|3.4j ) can be rewritten as follows, 



(3.9) 



(3.10) 



d d = (Ul) 2 aF>(m x a)(U R ) a2 , 



(3.11) 



where a = 1, 2 denotes the mass eigenstates as in Eq. 
is given by 



. The function D(m v a' 



g 2 TfbA 

(tan P + cot (3) — Im 



16n 



V m x a 



Ml 



K \ uiF U) v *L m « K . 



(3.12) 



ii 



Since we can express the function F(x) and the orthogonal matrices Ur and Ul in 
terms of eigenvalues m x a and M s , we obtain the EDM of Eq. ( |3.11|) 



as 



dd = ±- 



D{m x i) + D(m x 2) 
m Y i + m Y 2 



y/(m x i + m x 2) 2 - 2M^(1 =f sin 2(3) 



(3.13) 



i-k* I r,F>(m x i) - DimS 
+sgn(M 2 - \fi\)- 



m x 2 — m x i 



m x i -m x 2f - 2M^(1 ± sin 2,3) 
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where the upper (lower) sign corresponds to the positive (negative) sign of detM x . 

Since the off-diagonal elements of the A matrix are much smaller than the diag- 
onal ones in magnitude, let us first examine the effect other than the off-diagonal 
elements of the A matrix. Namely we tentatively assume that A is just a number 
A without the off-diagonal elements. Then the scalar-quark mass matrices have 
off-diagonal elements in generation indices only in the LL blocks. In this case, many 
diagrams become real and do not contribute to EDM, since the phases coming 
from the Cabibbo-Kobayashi-Maskawa matrix are canceled due to the hermiticity 
of mass matrices. Without the off-diagonal elements, we find the contribution of 
this chargino diagram to give the EDM of the order of 10 -31 ~ 10~ 33 e-cm, assuming 
SUSY mass parameters are of the order of 100 GeV. 

Next let us examine the effect of the off-diagonal elements of the matrix A. The 
scalar-quark mass matrices at low energies after the renormalization group running 
can be represented as Eq. ( |2.24j ) in the super KM basis. The off-diagonal elements 



of the A matrices are small in magnitude but have complex phases of 0(1) as in 
Eqs. ( |A.38|) and ( |A.41| ). Furthermore they provide new sources of flavor changing 



neutral current. These flavor changing currents spoil the cancellation of the KM 
phases in the one-loop diagrams for the EDM. Therefore the structure of the LL 
blocks of scalar-quark mass matrices are not important in this case. 

The presence of the off-diagonal elements of A generally helps to give a larger 
contribution to EDM. As shown in Fig. |3|, the contribution of the chargino diagram 
to the EDM in fact increases and becomes of the order of 10 -27 ~ 10~ 29 e • cm. In 
Fig. |3], we have chosen tan (3 = 10 and m x i = 50 GeV. One should note that EDM 
is approximately proportional to tan/5 for tan/? > 5. In order to see the allowed 
parameter region, we have plotted in Fig. f| the region of positive squared masses for 
scalar-quarks. The allowed parameter regions differ depending on M 2 > (Fig. 4a) 
or M 2 < (Fig. 4b). 

Since d-type-scalar-quarks instead of u-type-scalar-quarks are involved in the 
loop diagram, we find that the EDM d u of u-quarks is smaller than that of d-quarks 
and is of the order of 10~ 33 ~ 10~ 35 e-cm. 

Next let us examine the gluino contributions to the neutron EDM which is shown 
in Fig. |, 

^=^~^(g(^\ ) , (3.14) 



167r 2y3 3M 3 V W/rl/h 
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where and M3 are the strong interaction coupling constant and the mass of the 
gluino respectively and 



G{x) = — Vi {\ -x 2 + 2x\ogx). (3.15) 

o 1 J. X J 



Diagonalizing explicitly the d-type-scalar-quark mass matrix, we find that the 
EDM dd of the down quark have contributions from this gluino diagram of the order 
of 10~ 34 ~ 10~ 35 e-cm. The complex phases in off-diagonal elements of A terms are 
the major source of this contribution. 

In the SU(Q) quark model, the neutron EDM is given in terms of the single quark 
EDM as 

4 1 

d n = ^d d - -d u . (3.16) 

From the above results, we find that the neutron EDM from the single quark EDM 
is of the order of 10~ 27 ~ 10 -29 e-cm for tan (3 = 10. The chargino diagram of Fig. [I] 
is the dominant contribution, since A terms have off-diagonal elements. Although 
gluino diagrams can also contribute to EDM, its contribution is smaller than that of 
chargino, since complex phases are almost canceled if one takes the (1,1) component 
in generation indices. 

The contribution of chargino diagram in Fig. |l] was also examined in ref. [10 
with somewhat different values of parameters and a result of the order of 10 -30 e-cm 
for the EDM of neutron was reported. We have taken account of the evolution of 
gaugino masses in our renormalization group analysis. This may be a reason for the 
fact that we have obtained a larger value for the neutron EDM in comparison to ref. 

[Hoi. 



4. Effects of the transition electric dipole moment 



It has been shown that the neutron EDM is also induced in the nonsupersymmet- 
ric minimal standard model through the transition electric dipole moment (TEDM) 
of quarks within the baryon as illustrated in Fig. |5] |15|, [16|, [17]]. In this section, we 
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examine the TEDM, d — > s + 7, in the MSSM and calculate the neutron EDM based 
on their method. The matrix element of the electromagnetic current between 
single particle states of different Dirac particles can be written as follows, 



(Pf, s /b' M (°)|Pi> Si) = u f (p f ,s f ) 



F2(q 2 



+ l^q y 



F 3 (q 2 



rrif + rrii rrif + rrii 

+{q 2 Y - {m f - m^F^q 2 ) 



(4.1) 



+ 



q z 



rrif + rrii 



Y-qn l 5 F A (q 2 ) 



where rrii and rrif are initial and final Dirac fermion masses respectively. The domi- 
nant contribution to EDM comes from F 3 term as that of single quark case when the 
mass difference \rrif — rrii\ is small enough compared to the other relevant masses. 
In our particular case, the operator relevant to the TEDM can be written as 



K sd s(p s )ia^ v (p s -Pd) v lbd(p d ) 



(4.2) 



where 



Im k s 



IK 



sd. 



F 3 (0) 
m s + m d 



(4.3) 



is just the TEDM. 

The TEDM in the nonsupersymmetric minimal standard model has already been 
calculated under the assumption that m t <C M w (JTSl |IE|| . Recently it is more and 
more certain experimentally that the top quark mass is very large |12| . Therefore we 



must redo the calculation in the nonsupersymmetric case by taking into account of 
the large top quark mass. Assuming and m s to be small in comparison with the 
masses of the internal lines of the loop diagram, we have the following contribution 
to the TEDM of d -> s + 7, 



_sd 



G f 



V2(4vr) 



-K* ts K td (m d -m s ) f\ 



(4.4) 



where |E0| 



X 



6(x - l) 4 



;-8.x 3 + 3x 2 + 12x - 7 + 6x(3x - 2) logx). 



(4.5) 
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From the above equation, we obtain that the standard model gives 2 • 10~ 26 e-cm for 
the TEDM using the Cabibbo-Kobayashi-Maskawa matrix ( |A.35| ) - ( |A.37|) . Since 
the hadronic effects to convert the TEDM to the neutron EDM give a factor of 10 -7 



16], the contribution from the TEDM in the standard model to the neutron EDM 



becomes of the order of 10 33 e-cm. 

Next let us consider the TEDM of quarks in the MSSM. Similarly to the quark 
EDM, the chargino and gluino diagrams are most important. Since the bound state 
effects should be the same for super symmetric and nonsupersymmetric models, we 
shall calculate these diagrams for the quark TEDM and compare them with those 
of the standard model. 

In order to obtain the TEDM, we have only to replace di by sl in the analysis in 
the previous section. Similarly to the diagram in Fig. |1|, for instance, the chargino 
contribution is obtained from Eq. (|3.4|) as 



e 2 

Im K sd = ^(tan (3 + cot (3)—rrid Im 

I67H tr 




rn.K 



(4.6) 
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Let us note that the element we are interested in is not (2,1) but (1,2) for the LR 
block of the scalar-quark mass matrix, since the Higgsino coupling changes chirality 
such as u l R — > sl- By the same token, the gluino diagram for the quark TEDM 
corresponding to the diagram in Fig. § is obtained as 



Im k 



sd 



3M 3 V V 



Mq 2 



(4.7) 



RL/ 12 



In this case, we take the (1,2) element for the RL block because gauge couplings do 
not change chirality. 

We obtain the TEDM of the order of 10 -25 ~ 10~ 27 e-cm for the chargino con- 
tribution as shown in Fig. |6] and 10~ 28 ~ 10~ 30 e-cm for the gluino one. Thus we 
find that the TEDM in the MSSM is of the same order of magnitude as that in the 
nonsupersymmetric standard model. By combining our results with the hadronic 
matrix elements estimated already in the nonsupersymmetric case, we find that the 
resulting neutron EDM becomes of the order of 10 



-32 



10 



-31 



e-cm. 



Let us also consider the effects of other diagrams in the supersymmetric model. 
Besides the diagram which can be considered as the quark TEDM, there are diagrams 
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where two or more quarks within the neutron exchange the SUSY particles. The Im- 
parity conservation in the MSSM dictates that such diagrams must be box diagrams 
where all of the internal particles are SUSY ones. So the intermediate states at the 
hadronic level are of the order of mass scale of the SUSY particles which are much 
heavier than those in the diagrams that we have considered. Therefore the resulting 
neutron EDM is expected to be highly suppressed. 

Apart from the TEDM that we have considered in Fig. [5], another proposal for 
an effect involving many quarks was made in ref. ]17|. They considered the so- 
called penguin diagrams for the TEDM of quarks. They found that the conversion 
factor from TEDM to the neutron EDM is 10~ 4 instead of 10~ 7 due to a large 
QCD enhancement. Therefore this diagram gives the neutron EDM of the order 
of 10~ 30 e-cm which is the largest contribution in the nonsupersymmetric standard 
model. Since the enhancement due to the QCD corrections is the same order of 
magnitude in the MSSM as in the nonsupersymmetric standard model, we obtain 
that the neutron EDM is of the order of 10~ 29 ~ 10~ 31 e-cm. Moreover we find from 
Fig. ^ that the TEDM is of the same order of magnitude even for small values of A 
such as A < 1 TeV, whereas the single quark EDM becomes very small as shown in 
Fig. |3|. Therefore the penguin diagram with TEDM is more important than a single 
quark EDM for smaller values of A (A ~ 1 TeV). 
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work is supported in part by Grant-in- Aid for Scientific Research (T.I. and Y.M.) and 
(No. 05640334) (N.S.), and Grant-in-Aid for Scientific Research for Priority Areas 
(No. 06221222) (N.S.) from the Ministry of Education, Science and Culture. 



A Renormalization Group Equations 

We define the scaling variable t using the GUT scale Mqut and the relevant mo- 
mentum Q as 

t^log^p. (A.l) 

A tilde over couplings denotes a division by a factor An, namely a = a/(47r). The 
renormalization group equations for the gauge couplings and gaugino masses M, are 
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given as follows f21|j: 



a 



-= -ba 2 



Mi = -bt&iMi, 



where 



hi 



33 



bo = 1, fea 



and we denote derivative by t with a dot. The solutions are given as follows: 

6ti(t ) a G 



OCi(t) 



1 + bi(t - t )ai(to) 1 + bita G ' 



(A.2) 
(A.3) 

(A.4) 
(A.5) 



= 7—1-77 — 7 x - / , x = ^-rr\Mi[t Q ) = —oti(t), 
1 + bi(t - t )«i(*o) ou(ta) a G 



(A.6) 



where we assume the gauge coupling unification and the universal gaugino mass at 
the GUT scale: 



ai(0) = a 2 (0) = a 3 (0) = a G , 



(A.7) 



Mi(0) = M 2 (0) = M 3 (0) = M. 



(A. 



The renormalization group equations for Yukawa couplings are written as follows 



¥3-- 



2Y, 



2Y, 



2Y„ 



13 16 . 

— ai + 3a 2 + — 03 ) Y u 

7 „ 16 _ 

— CKi + 3a 2 + ya 3 ) F d 



3r u ylr u + 3Tr ( Y U Y[ ) Y u 



3Y d Y ] d Y d + 3Tr <Y d Y\ ) Y d 



9 _ \ ~ 

-«x + 3a 2 J Y e - 
3Tr (Y d Y[) Y e . 



3Y e YlY e + Tr ( Y e Y [ ) Y e 



(A.9) 



(A.10) 



(A.11) 



We now define hermitian matrices ctf (/ = u,d,e), 



& f = Y' f Y f . 



(A.12) 
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The renormalization group equations for d j are given as follows: 

• /13 Ifi \ 

2a u 



2a n 



a, 



f 13 16 \ 

2 f— 3i + 3a 2 + y a 3 J d u - 2 3d u 2 + 3Tr (d u ) a u 

-ct d dc u - ct u a d , (A. 13) 

(T 1 6 \ 

—ati + 3a 2 + Y^sj d d - 2 3d d 2 + 3Tr (d rf ) a d - 2Tr (d e ) a<i 

-«u«d - a<id u , 
9 



ai + 3a 2 a= e — 3a: e + Tr (d e ) a e — 3Tr (d d ) a e . 



(A.14) 
(A.15) 



The renormalization group equations for A's can be written as follows |2T|, [TD[ : 
2A„ = - 

5o;, u yl, u A u ol u qi^-A^ + A. u ol^ 2^4.^0:^ 



13 16 
-2 ( — aiMi + 3a 2 M 2 + y« 3 M 3 ) 1 - 2Tr (3A„d u ) 1 



(7 1 6 \ 

—aiMi + 3a 2 M 2 + — a 3 M 3 J 1 - 2Tr (A e d e + 3A d & d ) 1 

-Sd^A^ - A d d d - a u A d + A d d u - 2A u d u , 

A e = -2 f -fiiMi + 3a 2 M 2 ^j 1 - 2Tr (A e d e + 3A d d d ) 1 

OL P . 



(A.16) 
L 

(A.17) 
(A.18) 



Next we consider the renormalization group equations in the super KM basis, 
i.e., we rotate Yukawa couplings Y f by the same amount as M f in Eqs. (|2.20|) and 
( |2.21| ), and Af and ctf by the same amount as A f in Eq. (|2.23 ). In this basis the 
renormalization group equations are given as follows: 



2d, 



2d n 



a, 



2 ( — &i + 3<S 2 + ya 3 ) a u - 2 3d u 2 + 3Tr (d u ) d„ 



-KoL d K ] o. u - OL u Kot d K\ 

(7 1 6 \ 

— &i + 3a 2 + y a 3 J a d 



(A.19) 



3d/ + 3Tr (a d ) cx d - 2Tr (d e ) d d 



K ] a u K6L d - ct d K ] a u K } 

9 ~ \ r i 

-ai + 3a 2 ) cx. e — 3a e 2 + Tr (d e ) d e — 3Tr (d^) ct e 



(A.20) 
(A.21) 



2A U = -2 ( — ociMi + 3a 2 M 2 + y « 3 M 3 j 1 - 2Tr (3A M d M ) 1 

5d M A u - A u d u - Ka d K^A u + A u Ka d K^ - 2KA d a d K\ (A.22) 
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2A r , 



7 16 
-2 ( — a x M x + 3a 2 M 2 + —a 3 M 3 ) 1 - 2Tr (A e a e + 3A d a d ) 1 



5a d A d - A d a d - K^a u KA d + A d K^6t u K - 2K' A u a u K, (A.23) 
-2 ^^M! + 3a 2 M 2 ) 1 - 2Tr (A e a e + 3A d a d ) 1 



-5o; e -i4. e 



_A e Q! e . 



(A.24) 



In Eqs. ( |A.19| ) - ( |A.21| ), we give the initial conditions for ctf (f = u,d } e) at the 
scale of Z° boson mass M z which are obtained from the quark and lepton masses, 
i.e., 



a, 



OLd 



Oi, 




d ■ 



(A.25) 
(A.26) 
(A.27) 



at the scale Mz and we obtain solutions for ctf. The ctf are diagonal at the initial 
condition and then ct e is diagonal at all the scale but a u and ct d are not diagonal 
at the higher scale. 

In Eqs. ( |A.22p - ( |A.24p , we give the initial conditions for Af (/ = u,d, e) at the 
GUT scale M GUT which are universal, i.e., 



A u (0) = A d (0) = A e (0) = Al. 



(A.28) 



The Af are diagonal at the initial condition. Then A e is diagonal at all the scale, 
but A u and A d are not diagonal at the lower scale. It is convenient to separate Af 
into two parts. 

Af = A Lf + A Mf , (A.29) 



where f — u,d, e. The same renormalization group equations (|A.22|) - ( |A.24j ) as Af 
are valid for Auf which are proportional to the universal gaugino mass M . Their 
initial conditions at the GUT scale are 



A Mu (0) = A Md (0) = A Me (0) = 0. 



(A.30) 
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By deleting gaugino masses from Eqs. ( |A.22| ) - ( |A.24| ), we obtain linear equations 
for A Lf : 



2A 



La 



-2Tr (3A Lu a u ) 1 - 5dc u A Lu - A Lu ol u 
-K6t d K^A Lu + A Lu Ka. d K^ - 2KA Ld cx d K\ 
2A Ld = -2Tr {A Le a e + 3A Ld ct d ) 1 - 5a d A Ld - A Ld dc d 
-K^dt u KA Ld + A Ld K^a u K - 2K^A Lu a u K, 
A Le = -2Tr (A Le dt e + 3A Ld ct d ) 1 - 5a e A Le - A Le a. e . 

The initial conditions are the same as Af and are proportional to A, i.e., 



(A.31) 

(A.32) 
(A.33) 



A Lu (0) = A Ld (0) = A Le (0) = Al. 



(A.34) 



We consider rrit = 174 GeV [Q. The Cabibbo-Kobayashi-Maskawa matrix is 



parametrized as |22 



c 2 c 3 



c 2 s 3 



K 



-C1S3 - SiS2C 3 e i(5 C1C3 - s 1 s 2 s 3 e lS 



>e- t5 \ 
hc 2 



S1S3 - cis 2 c 3 e 



-S1C3 - cis 2 s 3 e 



<5 



s 2 e~ 

ClC 2 / 



(A.35) 



where 

Ci = cos^j, Sj = sinftj (i = 1,2,3). 
We use the following angles from phenomenological analyses |2"2" 



(A.36) 



9 1 = 0.043, fl 2 = 0.005 , 9 3 = 0.221, 5 = 0.86. 



(A.37) 



For tan/3 = 10 and A = 1 GeV, we obtain A^f at the scale Mz by solving the 
renormalization group equations: 



0.60 (-5.4 + il.8) ■ 10~ 6 (-3.6 + Z4.2) • 10" 5N 

A 1 „ = I (-5.4 - zl.8) ■ 10" 6 0.60 -4.8 ■ 10" 4 

-3.6-i4.2) ■ 10~ 5 -4.8 ■ 10" 4 0.19 



(A.38) 



A Ld = 



I 0.96 

(3.7-i2.2) ■ 10~ 5 
\ (-8.3 + i4.9) • 10- 4 



(3.7 + z2.2) 
0.96 
(5.7 + i0.1) 



10" 



10" 



(-8.3-14.9) • 10" 4 * 
(5.7-iO.l) • 10~ 3 
0.79 



(A.39) 
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'0.96 \ 
v L , i 0.96 

0.95/ 



(A.40) 



For tan (3 = 10 and M = 1 GeV, we obtain A M f at the scale M z : 



/ -2.9 

(l.l + i0.4) • 10" 5 
\(7.6 + i8.9) • 10~ 5 



(1.1 -i0.4) • 10- £ 
-2.9 
0.0010 



(7.6-i8.9) • 10" 
0.0010 
-2.1 



(A.41) 



/ -3.7 

(-8.0 + i4.7) • 10- 5 
V (1.8 + il.l) • 10- 3 



(-8.0-i4.7) • 10~ 5 
-3.7 
-0.012 



(1.8- 




(A.42) 



t-Me 



(- 




(A.43) 
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Figure captions 



Fig. 1 The chargino contribution to the quark EDM involving Higgsino couplings. 
The arrow on each line stands for the chirality of the particle and the cross 
for mass insertion. 

Fig. 2 The gluino contribution to the quark EDM. 

Fig. 3 For the case det M x > 0, EDM of the down quark is plotted as a function 
of the mass m x 2 of the heavier chargino for various values of A in the case 
(a) M 2 + // > and (b) M 2 + /i < 0. We have chosen parameters: universal 
scalar-mass m = 100 GeV, tan (3 = 10, and the mass of the light chargino 
m x i = 50 GeV. 

Fig. 4 The region of parameters for the squared mass of scalar-quarks to be positive, 
for (a) M 2 + fi > and (b) M 2 — < 0. The boundary of vanishing mass of 
u-type-scalar-quark is represented by solid lines, and that of d-type by dashed 
lines. In the case (a), allowed region for m = 100 GeV, /i < 0, and tan/3 = 10 
is denoted by shaded area. For tan/5 = 100, the allowed region is the right 
of dotted line bounded by two solid lines. In the case (b), allowed region for 
/i > is shown in the upper half plane and that for \x < in the lower half. 
Allowed region for m = 100 GeV and tan (3 = 10 is denoted by shaded area. 

Fig. 5 Contribution to the neutron EDM through the TEDM. The diagrams like 
Fig. [I] or Fig. |2] is inserted in the blob. 

Fig. 6 For the case det M x > 0, TEDM (d — > s-y) is plotted as a function of the mass 
m x 2 of the heavier chargino for various values of A in the case (a) M 2 + fi > 
and (b) M 2 + fi < 0. We have chosen parameters: universal scalar mass 
m = 100 GeV, tan (3 = 10, and the mass of the light chargino m x i = 50 GeV. 
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4^ 
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Figure 1: The chargino contribution to the quark EDM involving Higgsino cou- 
plings. The arrow on each line stands for the chirality of the particle and the cross 
for mass insertion. 




Figure 2: The gluino contribution to the quark EDM. 
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Figure 3: For the case det M x > 0, EDM of the down quark is plotted as a function 
of the mass m x 2 of the heavier chargino for various values of A in the case (a) 
M 2 + pi > and (b) M 2 + pi < 0. We have chosen parameters: universal scalar mass 
m = 100 GeV, tan/3 = 10, and the mass of the light chargino m x i = 50 GeV. 
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Figure 4: The region of parameters for the squared mass of scalar-quarks to be 
positive, for (a) M 2 + /i > and (b) M 2 — \fx\ < 0. The boundary of vanishing mass 
of u-type-scalar-quark is represented by solid lines, and that of d-type by dashed 
lines. In the case (a), allowed region for m = 100 GeV, fi < 0, and tan/5 = 10 is 
denoted by shaded area. For tan/? = 100, the allowed region is the right of dotted 
line bounded by two solid lines. In the case (b), allowed region for \x > is shown 
in the upper half plane and that for ji < in the lower half. Allowed region for 
m = 100 GeV and tan f3 — 10 is denoted by shaded area. 
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Figure 6: For the case det M x > 0, TEDM (d — > sj) is plotted as a function of the 
mass m x 2 of the heavier chargino for various values of A in the case (a) M 2 + \i > 
and (b) M 2 + \x < 0. We have chosen parameters: universal scalar mass m = 100 
GeV, tan/3 = 10, and the mass of the light chargino m x \ = 50 GeV. 
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